Search for: All records

ABSTRACT Voltage‐dependent anion channel (VDAC) is the primary conduit for regulated passage of ions and metabolites into and out of a mitochondrion. Calculating the solvation free energy for VDAC is crucial for understanding its stability, function, and interactions within the cellular environment. In this article, numerical schemes for computing the total solvation free energy for VDAC—comprising electrostatic, ideal gas, and excess free energies plus the nonpolar energy—are developed based on a nonuniform size modified Poisson–Boltzmann ion channel (nuSMPBIC) finite element solver along with tetrahedral meshes for VDAC proteins. The current mesh generation package is also updated to improve mesh quality and accelerate mesh generation. A VDAC Solvation Free Energy Calculation (VSFEC) package is then created by integrating these schemes with the updated mesh package, the nuSMPBIC finite element package, the PDB2PQR package, and the OPM database, as well as one uniform SMPBIC finite element package and one Poisson–Boltzmann ion channel (PBIC) finite element package. With the VSFEC package, many numerical experiments are made using six VDAC proteins, eight ionic solutions containing up to four ionic species, including ATP⁴⁻and Ca²⁺, two reference states, different boundary values, and different permittivity constants. The test results underscore the importance of considering nonuniform ionic size effects to explore the varying patterns of the total solvation free energy, and demonstrate the high performance of the VSFEC package for VDAC solvation free energy calculation. 

Abstract Integrated computational modeling provides a mechanistic and quantitative framework to characterize alterations in mitochondrial respiration and bioenergetics in response to different metabolic substrates in-silico. These alterations play critical roles in the pathogenesis of diseases affecting metabolically active organs such as heart and kidney. Therefore, the present study aimed to develop and validate thermodynamically constrained integrated computational models of mitochondrial respiration and bioenergetics in the heart and kidney cortex and outer medulla (OM). The models incorporated the kinetics of major biochemical reactions and transport processes as well as regulatory mechanisms in the mitochondria of these tissues. Intrinsic model parameters such as Michaelis–Menten constants were fixed at previously estimated values, while extrinsic model parameters such as maximal reaction and transport velocities were estimated separately for each tissue. This was achieved by fitting the model solutions to our recently published respirometry data measured in isolated rat heart and kidney cortex and OM mitochondria utilizing various NADH- and FADH2-linked metabolic substrates. The models were validated by predicting additional respirometry and bioenergetics data, which were not used for estimating the extrinsic model parameters. The models were able to predict tissue-specific and substrate-dependent mitochondrial emergent metabolic system properties such as redox states, enzyme and transporter fluxes, metabolite concentrations, membrane potential, and respiratory control index under diverse physiological and pathological conditions. The models were also able to quantitatively characterize differential regulations of NADH- and FADH2-linked metabolic pathways, which contribute differently toward regulations of oxidative phosphorylation and ATP synthesis in the heart and kidney cortex and OM mitochondria.